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Abstract 

We study the probability of two Brownian particles to meet before one of them exits 
I a finite interval. We obtain an explicit expression for the probability as a function of 

^ ' the initial distance of the two particles using the Weierstrass elliptic function. We also 

CSl ■ find the law of the meeting location. Brownian simulations show the accuracy of our 

P5 . analysis. Finally, we discuss some applications to the probability that a double strand 

^ ' DNA break repairs in confined environments. 



0^ 



0^ 
O 



keywords: Brownian motion, conformal mapping, Weierstrass function. 



Introduction 

J3 ' The problem of coalescence and clustering in an open space has been considered by Chan- 

drasekar P (see also [6] for a review). Little work has been dedicated to the case of a 
competition between the coalescence of Brownian independent particles and the possible es- 
cape at the boundary of a domain where the particles are absorbed. This situation is however 
reminiscent of many biophysical problems. For example, the probability of a correct repair of 
' a broken DNA molecule inside the nucleus. One mode of repair is known as non-homologous 

end-joining [5], which depends crucially on the initial distance between the two free DNA 
^ ! strands: either the branches meet or one of them can curl up before and then the probability 

O I to connect is very low (almost zero). Here we consider the drastic simplification that the 

motion of the DNA molecule tip can be approximated as a one dimensional Brownian motion 
(see the discussion). 

^ ! We consider the motion of two independent Brownian particles Xi(t), X2{t) inside an 

^ I interval [a,b], (a < b) with the following rules: when the two particles meet, they coalesce 

into a single one subjected to a Brownian motion. We compute the probability Pm that 
the two particles meet before one of them hits the boundary of the interval and obtain an 
explicit expression for the probability of the two particles clustering, as a function of the 
initial position. When the initial points are a < xi < 0:2 < we obtain that 

Pm{x,,x,) = ^3mlog^ i^^^T^) ' 
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where 3m is the imaginary part, ^ is the Weierstrass elhptic function defined by equation 
dHD, L = b — a, Z = X2 + y/'^Xi, and 

dx 

J = 5.244115106 (2) 

[x{x — 1)] * 

is a universal number defined by an elhptic integral. We further obtain the probability 
distribution of their meeting point. Finally, the analytical formulas are compared with 
Brownian simulations, where we gain the information about the variance. The role of the 
Weierstrass elliptic function is quite surprising here and really comes from the method of 
conformal mapping. We wonder if our result can be recovered from elementary probability 
arguments. A Brownian interpretation of an elliptic integral was given in [3]. 




Formulation 

The dynamics of each particle is given for z = 1, 2 by 

dXi = y^2lTfdwi (3) 

where Df is the diffusion constant and Wi,W2 are two Brownian motions of unit variance. 
We are interested in the probability Pm that the two particles meet before one of them exits 
the interval [a, b]. If we consider the two random times 

Ti = inf{t > 0, Xi{t) = a or X2{t) = b, ^i(O) = Xi and X2{0) = X2, xi < X2} 
T2 = inf{t > 0, Xi{t) = X2{t),Xi{0) = xi and X2(0) = X2, xi < X2}, 

then for x = {xi, X2), the probability 

PMix) = Pr{T2 < n\x} (4) 

satisfies the Laplace equation 

APMix) = OioTXET (5) 
Pm{x) = Hot X e D 
Pm{x) = for xedT-D 

where T is a right-angled triangle with vertices a,b,b + a\^—l. D is the side joining a to 
b + a\/^. Indeed, 

Pm{x) = Pr{X{T) = ye T\X{0) = x) = I G{x, y)dSy (6) 

Jt 

where r is the first exit time and the Green function G is solution of (see [^) 

AG{x,y) = OfoTxeT (7) 
G{x,y) = 6y{x) for X e D 
G{x,y) = ioT X e dT - D. 
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G{x, y) is the probability density function to exit in y E D when the particle starts initially 
in X (see also ch. 15, p. 192 of reference [2] for another proof). We shall derive an explicit 
expression of the encounter probability P. To solve equation, we shall use the invert of a 
Schwarz-Christoffel mapping to map the triangle into the upper complex half-plan H. By 
using the explicit the solution of the Laplace equation in H, we will find the solution of ([5]). 
It turns out that the Schwarz-Christoffel mapping of interest is a Weierstrass function. 



Analytical derivation of the encounter probability 

It will be convenient to do all our computations for the choice a = 0,b = u, where uj is 
defined in [2l Note that u and uy^—l are a pair of fundamental periods for the Weierstrass 
^ function with parameters (72 = 1 and = 0, [4j defined by 

= 4^33 - g,^ - g,, (8) 

a function we are going to use in the following. It is a matter of a dilation to deduce the 
results for the case of general a, b from our special case. 

In the spaces H = {z E Cl'^mz > 0} and H = {z E C\'^mz > 0}, we consider f : H ^ <C 

3 

to be the branch of [z{z — 1)] * defined as follows: for z E H, we set 



where Vq = \z\, ri = \z - 

3 

f{z) = (rori)4 exp 
f has a continuous extension to H : 



Tq exp (^^oV— 1) and z — 1 = ri exp (^i V— l) 
- 1\, < 9q,9i < TT. In that case, 

^0 + ^1' 



1)] 



[X[X 

^[\x{x - 1)1]^ 



[x{l — x)] 



if X > 1 
if X < 
if < X < 1 



We shall now define F : H ^ C hj: 



(9) 



(10) 



:ii^ 



dz 

W) 



:i2) 



The Schwarz'reflection lemma shows that F is a conformal mapping of H onto the triangle 



+00 



T in C having as vertices 0,a;, (l + a/— l) uj where uj = f f\x)dx (see figure [T]). F maps 

1 

1, 00, onto 0, UJ, (1 + v^— T) UJ respectively and the half-line [1, +00] onto the segment [0, uj] 

of the real axis, the half line [—00, 0] onto {uj + tuj\/—l\Q < ^ < 1} and the segment [0, 1] of 
the real axis onto {(l + \/— l) (1 — t)u\Q < t < 1}- Moreover F is conformal on H. 
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To compute the function F given by (fT2|) . or more precisely its inverse we introduce the 
following transformation: z E H ^ p = ifi{z) G C, 



where the square root is the one such that > 0. With the above notations, p = 
ly^exp [^ ^o-gi ij , tp is a homeomorphism of H onto the quadrant Q = {p E C\^ep > 
0, '^mp < 0} and (p is conformal on Q = {p E C\^ep > 0, Qp < 0}. On the space H, 

f{z)dz = ^*{g{p)dp) (14) 

where 

9{P) = 1^ „ 

^J^P^ - p 



and a/ — p is the branch in Q which is real positive on ]i,+oo[. 99 maps [l,+oo] onto 
[+00, i], [0,1] onto -a/^[0,+cx)] and [-00, 0] onto [|,0]. 

Relation ([H]) in ([12]) implies that: 



/»+00 /"H-OO -1 

i^(C) = / ^(P)^ip = / r—, dp. (15) 

We recall that the Weierstrass elliptic function [4J,^ with (72 = 1, 5*3 = is defined as the 
C = ^(0 such that 



Relation f[T5|) implies that: 



Equivalently, ([T7|) implies that: 



? = / dp. (16) 



^(0 = ^(^) (17) 



Hence we have for Z E T : 



F-\Z) 



To compute ([5]), we shall find a harmonic function m with the boundary condition 



u\]id,uj\VJ {uj + {tuj^/^\Q < t<l} = 

M|{(l + V^)tw| -1 < t<0} = l. 



This harmonic function can be expressed using the Harmonic function v in H with the 
boundary conditions 

on [-cx),0[U]l,+oo] 

(18) 

1 on [0,1] 



It is given for Z G T by: 



u{Z) = v{F-\Z)) = V 



Actually, the function v is given for z & H : hj 



4<P( ^ 



V8' 



v{z) = —53m log 



z-1 



TT 



where log - — is the branch on H —{0, 1} that is real for z g]1, +oo[ i.e.: 



log^ = log^ + v^(^i-^o) 
z ri 



Finally, we obtain for Z eT the expression 
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u{Z) = — $5mlogq3(— ) 



TT 



(19) 



(20) 



(21) 



where log is the branch on H —{0} which is real on ]0,+oo[. This result can be seen as a 
generalization of Spitzer's law concerning the winding of a two dimensional Brownian motion 

m. 



Comparison with Simulations 

When the initial positions of the Brownian particles are < xi < X2 < 1, the scaled meeting 
probability is given by 

Pm{x„x,) = zlQmlog^{u{^^±^^^)), (22) 

TT V8 

where u is defined in equation [2l In figure [2], we present the graph of the probability Pm of 
forming a cluster, plotted as a function of the initial positions Xi with the restriction that 
Xi < X2- We present two simulations: in the first one, we fix the point X2 at the middle 
of the interval (x2 = 0.5) and the other graph is obtained for a point X2 chosen very close 
to the boundary (x2 = 0.99). As can be observed, the shape of the encounter probability 
changes drastically. We have superimposed in figure [21 the Brownian simulations (mean and 
variance) with the analytical solution. 

We remark that the probability to meet does not depend on the diffusion constant. 
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(1+i) 0) 




" 10 

Figure 1: Transformation F from the triangle to the upper complex plane. We position the 
boundary condition for the Laplace equation on the associated part of the boundary 




0.2 0.4 0.6 0.8 0.99 0.1 0.2 0.3 OA 0.5 0.6 

Xl "1 

Figure 2: Probability Pm for two Brownian particles to meet before one of them escapes the 
interval [0,1]. We compare the analytical solution Eq. [5] with Brownian simulations. For each 
position, we averaged 2000 realizations. The variance is presented as error bars. On the right, we 
start at a middle point of the interval (x2 = 0.5) and plot the probability as a function of xi (with 
Xl < X2 for the second point, while on the left, we chose for X2 a point very close to the boundary 
(x2 = 0.99). The effect of the boundary layer appears clearly for {x2 = 0.99). 
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The position of encounter 



To satisfy our last curiosity, we finish with the computation of the probabihty p{Z] E), Z = 
X2 + ^/—lXl for the two particles, starting at positions {xi,X2), < Xi,X2 < uj" to coalesce 
in a measurable subset E of [0, a;]. Given E e [0,a;], the function Z eT ^ p{Z; E) e [0, 1] 
is harmonic in the interior of T and 

mE)-{\ (23) 

Using the conformal transformation T ^ H, 

Z^z^ i^iZ) = / ^ / , (24) 

it is sufficient to compute for any measurable subset M of [0,1], the function: z e H ^ 
P{z; M) e [0, 1] having the following properties: 

1. P is harmonic in H. 
2. 

3. for any z e H, M e mt([0, 1]) ^P{z;M) is a measure on the a-algebra 9Jl([0, 1]) of 
measurable subsets of [0, 1]. Then, 

p{Z;A)^P{^{Z),i;{A)). 

We shall remark that condition 3) above shows that to determine P, it is sufficient to 
compute P{z, .) for M a closed interval [a,P] of [0, 1]. This is similar to the determination 
of the function v above. Hence, 

P(z;[a,l3]) = -^mlog^^. (26) 
TT z — a 

This expression shows that ior z E H, M e Tl{[Q, 1]) -^P{z; M) has a density D{z] a),z E 
H^a E [0, 1], with respect to the Lebesgue measure on [0, 1] : 

d P 1 X 
D{z\a) = -^{z\ [a,(3\) = Q=m- 



da TT z — a 



Hence p{Z; E) has a density with respect to the arc length on the segment { (l + \/— T) (1 — ^) <^|0 < 
t < 1}, which, for A e {(1 + \/^ (1 -t)uj\Q<t< 1}, is equal to : 



d{Z-A)^D{^{Z),i^{A)) 



(27) 
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Recall that ip maps the segment {(l + -\/— T) (1 — t)ci;|0<t<l} diffeomorphically onto the 
interval [0, 1]. 



2 ' 



(28) 



where 



and using the different equation satisfied by we obtain that 



(29) 



dijj 



^( 



(30) 



Finally the meeting density function (!27j) is given explicitly by 
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d(Z; A) = — 

TT 



/ / z_ 



(31) 



where Z & T and ^4 G {^^(1 — t){l + V— 1)|0 < t < 1}. If £" is a measurable subset of the 
state space [0,uj] then the probability for two particles starting at position Z to meet in E 
is given by 



p{Z; E) = V2 d{Z; (a + aV^))dc 
Je 



(32) 



which conclude this part. To finish, we extend our formula for an initial segment [a,b]. In 
that case, the probability to meet before escape scales into: 



P{Z) = — 9mlog<p 



TT 



uj{Z — a) 



(33) 



where Z E T' the triangle in C with vertices a,b,b + {b — a) V— 1. In the limit of L large, 
using the double pole expansion of ^ at zero, for Z in a neighborhood of 0, we have 



q3(Z) = ^ + 0(Z), 
thus for large L, for Z = '^^^^^ ~ equation ( l33l) becomes 

P{Z) ^L^oo —QmlogZ = — arctan ( — 

TT TT \X2 



(34) 



(35) 



We shall remark that this law is twice the one of a Cauchy variable at time 1, that is 
Prob{\Ci\ < f^}- This suggests that this asymptotic result might be recovered by elementary 



considerations on the Brownian motion. Finally, the meeting probability density function at 
the point A is now given by: 



d{Z;A) = 




To finish, wc shall provide the asymptotic for d{Z; A) for large L. Using the meromorphic 
property of ^, we obtain the following expansion of its derivative at the origin: 

^'(Z) = -|- + 0(1), (36) 



Thus, 



oj{A-a) \ 

lVs J 



2LVS 



(37) 



<p/ f i>j{Z-a) \ ^, f uj{A - a) 



4(W8)^ 
(u^AZ)^ 



(38) 



^ 'a^(Z-a)\ (^^fu{A-a) 



LV8\ f LVs' 
uZ I I ujA 



(39) 



u{Z — a) 



J \ujA J \ujZ 



(40) 



d{Z;A) 



ujA 



ttL 
16\A\^ 

71 

16\A\^ 

TT 



2LV8 



4(L 



12 



{AZ) 



-6 



9m 



{AZ) 



-6 



1 



Z4 -^4 



(41) 



The mean conditional time for a collision before exit 

We shall continue here with the expression for the mean conditional time Tm{x) to meet 
before one of the particles escape. The conditional time Tm^x) to hit the diagonal of the 
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triangle before the other sides is associated with the conditional process X* solution of the 
stochastic differential equation [2], 

dX*it) = 2D^^$^^-I^dt + V2DdW, 
p{X*{t)) 

where p is the probability (!33|) . Tm{x) satisfies Dynkin's equation [7] 

DpATm + 2DVTm " = —p in T, 

= OonD, (42) 

where D is the diagonal (there are no conditions on the other side). Thus w = pr^ satisfies: 

DAw = —p in T, 

w = Oon dT, (43) 

We recall that the solution of the Dirichlet problem 

Au = k in H, 
M = on dH, (44) 

is 

u{y) = / G{y,zi)k{zi)dzi (45) 
Jh 

where the Green function G is given by 

G{y,z^) = ^\JJi^. (46) 
2tt \y-zi\ 

Using now that for any conformal transformation 0, 

Aiuo<P){z) = \^'{z)\'Au{<Piz)) = -\<P'iz)\M<P{^)), (47) 

we obtain 

{uo<P)iz) = -^J^ \<P'iz,)\M<P{z^))Giz - z,)dz, (48) 
with Z = (f){z) = F{z), 

u{Z) = -^j^\F'{z^)\'p{F{z,))G{F-\Z),z,)dz, (49) 



l/jFf (F-(ZO)p^)G(F-(Z),F-(ZO)^;J^ 



(50) 
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Discussion 



The dynamics of double strand DNA (dsDNA) break is a fundamental step of the repair 
process. There are no direct experimental measurements yet of the dynamics of dsDNA 
ends (telomere) in the confined nucleus environment, thus giving a fundamental role of 
the theory in understanding the physics of motion, leaving aside the molecular machinery 
involved in the repair process. Since the general picture of telomere motions is not clear, we 
have presented here a very coarse analysis based on Brownian motion, which can be seen 
however as a drastic simplification of polymer motion in a confined environment, restricted 
by the nuclear crowding, including histones, the remaining DNA organization, nucleoli and 
many other nuclear components. When the microdomain surrounding the dsDNA break is 
sufficiently narrow (a long strip of length 1), using the Rouse model for the polymer, with a 
persistence length Iq, we can distinguish two cases: I < Iq or I > Iq. The polymer is modeled 
as an ordered string of beads, each being connected to its next neighbors by a spring of 
elasticity constant k. The mean length between the beads is Iq. The motion of a string is 
governed by a multi-dimensional Langevin equation, the potential of which is due to the 
elastic forces. For a bead at position Xj, the motion is described by the Smoluchowski limit 
of the Langevin equation (i=l..N) 

Xi + VU{xi,Xi+i,Xi_i) = V2Dwi (51) 

where D is the diffusion constant, Wi are 5— correlated Brownian motion of variance 1, and 
the potential U is 

{U{xi,Xi-i) + U{xi,Xi+i) if i=2..N-l, 

(52) 
U {xj,Xj+i) iij = lorj = N-l 

and 

U (xi, Xi-^-l) = fc( — — — ^ol^^i+i ~ 3;j|) (53) 

When I < Iq, the polymer cannot collapse and the DNA ends may be approximated by 
the one dimensional motion of the polymer chain. This approximation is so restrictive that 
the polymer relaxes to its equilibrium and the two ends meet with probability one. When 
I > Iq, there are two final possibilities: starting at an initial position, either the two branches 
touch or curl up and then they will not be able to be repaired in a reasonable time. We 
have restricted our analysis to an one- dimensional Brownian motion. A full analysis of 
this phenomenon is difficult and we shall discuss now some ideas to address it. First our 
analysis is relevant for short DNA fragments between two neighboring nucleosomes where 
we assimilate the break location to the center of mass of the polymer, whose motion is 
Brownian. However, the computation we presented here of the probability to bind before 
escape cannot be generalized easily to dimensions 2 or 3 because it depends heavily on 
conformal mappings. To generalize our result, it is possible to use the Rouse model for a 
polymer and estimate the probability that the two ends of the dsDNA break meet for the first 
time before one of them collapses. Similarly, estimating the probability that the two ends 
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meet before a given time would also be relevant. These questions are much more difficult 
to address compared to our analysis. However, a first step using simulations would be to 
estimate the mean ffist passage time to one of the polymer end to reach a small hole. This 
is already a nontrivial generalization of the small hole theory, because the small hole is not 
small anymore. We hope that our analysis will help to understand better the mechanisms 
of repair processes occurring in the extremophilic bacterium Radiodurans, where radiations 
are known to produce a nuclear phase transition, leading to a restriction of the space and 
thus increasing the probability of DNA repair [5J . 
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